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Devices based on ultracold atoms moving in an accelerating optical lattice or double-well potential 
are a promising tool for precise measurements of fundamental physical constants as well as for the 
construction of sensors. Here, we carefully analyze the model of a couple of BECs separated by a 
barrier in an accelerated field and we show how the observable quantities, mainly the period of the 
beating motion or of the phase-shift, are related to the physical parameters of the model as well as 
to the energy of the initial state. 

PACs numbers: 03.75.Kk, 37.10.Jk, 03.75.Dg 


Laser-cooled atoms have drawn a lot of attention as 
for potential applications to interferometry and high- 
precision measurements, from the determination of grav¬ 
itational constants to geophysical applications , see 
also [51 |B] for a recent review. The idea of using cold 
atoms moving in an accelerating optical lattice izHn] has 
open the field to multiple applications. For instance, by 
means of method proposed by Glade et al [12] a pre¬ 
cise measurement of the Earth’s gravitational accelera¬ 
tion constant g were performed [H [H; the obtained 
results had a very high precision and only a tiny discrep¬ 
ancy between g measured by a Raman interferometry on 
laser-cooled atoms and a classical gravimeter resulted, in 
fact the absolute relative uncertainty f^g/g turns out to 
be of order 3 x 10“®. 

More recently, a value for the constant g has been 
measured using ultracold strontium atoms confined in 
an amplitude-modulated vertical optical lattice m , im¬ 
proving a previous result m by using a larger number of 
atoms and reducing the initial temperature of the sample. 
Determination of g has been obtained by measuring the 
frequency vb of the Bloch oscillations of the atoms in the 
vertical optical lattice and recalling that vb = mgd/2TTfi, 
where m is the mass of the Strontium atom, h is the 
Planck constant and d is the lattice period. Since Bloch 
oscillations only occur for an one-body particle in a pe¬ 
riodic field and under the effect of a Stark potential then 
has been chosen, in the experiment above, a particular 
Strontium’s isotope ®®S'r; in fact the scattering length 
as of atoms ®®S'r is very small and thus it can be as¬ 
sumed that the effects of the atomic binary interactions 
are negligible. The obtained value for the constant g 
was consistent with the previous one but was affected by 
a larger relative uncertainty of order 6 x 10“®, because 
of a larger scattering in repeated measurements, mainly 
due to the initial position instability of the trap. Such a 
technique is also proposed to measure surface forces m- 

On the other side, new technologies enable the con¬ 
struction of simple coherent matter-wave beam splitter 
based on atom chips. These devices have been shown to 
be capable of trapping and guiding ultracold atoms on a 


microscale; BECs can be efficiently created in such small 
devices and coherent quantum phenomena have been ob¬ 
served. Interferometers based on a microchip can be 
widely used as highly sensitive devices because they al¬ 
low measurement of quantum phases. Technologically, 
chip-based atom interferometers promise to be very use¬ 
ful as inertial and gravitational field sensor provided that 
the quantum evolution of the matter waves is not per¬ 
turbed by the splitting process. It has been seen that 
in such a device a BEG cloud up to 10^ Rubidium-87 
atoms can be split in two clouds inducing a double-well 
trapping potential phase-preserving m- By means of 
such a devices a measurement of the Earth’s acceleration 
constant g has been performed with relative uncertainty 
of order 2 x 10“'* [T5j. It is well known that one of the 
most relevant physical effect in a double-well model is the 
so called heating motion between the two wells; hence, in 
principle one can measure the beating period of the BEG 
in an accelerated double-well potential and then obtain 
the value of the gravitational constant, as done for BEGs 
in an accelerated optical lattice. However, we would 
remark that in the case of isotopes the scattering 
length Os is not small and thus binary interactions must 
be effectively taken into account if one want to relate the 
Earth’s gravitational constant g with the beating motion 
of the two BEG’s clouds between the two wells. There¬ 
fore, in order to improve the analysis of the experimental 
output it is necessary to have a more complete under¬ 
standing of the underlying theory of two BEGs separate 
by a asymmetrical barrier. 

The aim of this paper is to provide a solid theoretical 
ground for a BEGs in a double well potentials under the 
effect of the gravity force, where an explicit formula con¬ 
necting the physical parameters, and in particular the 
Earth’s acceleration constant g, with the period of the 
observed beating motion between the two wells and of 
the difference of phase of two condensates. By means of 
such a result we expect that the relative uncertainty of 
the experimental results obtained for BEGs in chips may 
be improved. Indeed, in such a framework the mea¬ 
sure of the period of the difference of the phases between 
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the two condensates gives a precise value for the Earth’s 
acceleration constant g. We would underline that our 
analysis will be useful even as a model for a.c. Josephson 
effects in BECs m 

Here, we consider a simple model of BEC trapped in 
a double-well potential under the effect of a Stark po¬ 
tential, the dynamics along the direction of the grav¬ 
ity force is described by the one-dimensional Gross- 
Pitaevskii equation (GPE) 

I V’(a;,0) = ^o(4 

where V is the double-well trapping potential V, v = mg 
is the strength of the Stark potential and the nonlinearity 
is given by e = , with N the total atom number. 

Os is the scattering length, g is the gravity acceleration 
and m is the atom mass. The BEC wavefunction tp is 
normalized to one. 

By assuming the two-level approximation [niisa then 
the normalized BEC wave function tp can be written as 

ipix,t) = [aR{t)(pR{x) + aLit)(pL{x)] 

where aR^Rpt) are two complex valued functions depend¬ 
ing on the time t satisfying 

|afl(i)P + |aL(t)P = l; 

the vector (pR (resp. (Pr) corresponds to the ground state 
of the corresponding isolated right hand side (resp. left 
hand side) trap with associated energy H. 

It is well known that the solution to the unperturbed 
problem Q , where e = 0 and v = 0 and when the state is 
initially prepared on the first two ground states, exhibits 
a beating motion with period ^ independent of the ini¬ 
tial wave function "tpo, where uj = is half of the 

the splitting between the two onsite energies E± = 

Hence, ^ plays the role of unit of time and it is natural 
to introduce the (adimensional) slow time 



The amplitudes OR^Rpr) obey the nonlinear two-mode 
dynamical system given by (hereafter ’ = -^) 

f = -«z, + + POfl 

[iaf = -aR + f]\aR\‘^aR - paR ’ ^ 

where g and p are the adimensional quantities defined as 

l•+oo 

V=- \‘PRix)\'^dx 

^ J —oo 

and 

p = - x\pR{x)\'^dx. 

^ J — C30 

The unperturbed solution (for 77 = p = 0) of the two-level 
approximation has periodic solution with period T = tt. 


If we set ttR^Rir) = where qr^r S [0,1] 

are such that 9 fl(T)-|-g|^(T) = I, then the previous system 
([^ takes the Hamiltonian form 


O' = 

^ ~ do 


(4) 


where 0 '■= Or — 9r is the phase shift and z := qlp — q\ 
is the imbalance function between the two condensates, 
with Hamiltonian function 

"H = —2\/l — cos 0 + ^ 77(1 -p z^) -P 2pz . (5) 

We should remark that Q is invariant with respect to 
the change of the sign of e and p; more precisely, if p < 0 
then we can switch to the case of p > 0 by changing the 
signs z —>■ —z and 0 —>■ —9. Similarly, if 77 < 0 then we 
can switch to the case 77 > 0 by z —?> —z and 9 ^ 9 + n. 

It is a remarkable fact that equation admits explicit 
periodical solutions, and that the period T of the imbal¬ 
ance function z(t), as well as of the phase shift 9{t), 
can be explicitly computed as function of the parame¬ 
ters ?7 and p as well as of the initial wave function [23) . 
Therefore, in principle, if one experimentally measure the 
inversion frequency then one can obtain a precise value 
for the acceleration constant g. In fact, let us denote by 
E the energy value of the Hamiltonian T-L on the initial 
state: E := 'H{zq,9q). Then, from @ and (H we have 
that z(t) is a solution to the following ordinary differen¬ 
tial equation of first order: 


{z'Y = az^ + 9z^ + cz'^ + dz + e ( 6 ) 


where we set a = —\g^, b = —2gp, c = Eg —A— ^ —4p^, 
d = 4Ep — 2gp and e = Eg — E'^ +A — Equation ^ 
has solution given by means of the Weierstrass’s elliptic 
function 'P{t\ g 2 , gs) with parameters 


92 ■■= ae- -bd+ —c 


1 .2 


1 1 ..2 1 ... 1 3 

-eac — —ad + —dbc — c 
6 16 48 216 


The Weierstrass’s elliptic function 7 ^(t; 52 , 53 ) is a dou¬ 
bly periodic function which real period coincides with 
the period T of the phase shift and of the imbalance 
functions. In order to compute the real period T let e^, 
j = 1, 2, 3, be the roots of the trinomial 4s^— 525 — 53 ; and 
let 6 = pf — 275 I. If (5 > 0 then S i?, 63 < 62 < 0 < Ci, 
and 


2K{k) ^ _ 62 - 63 

\/ei — 63 ’ 61 — 63 


where K denotes the complete elliptic integral defined as 


K{k)= f [(1 - s2)(l - fcs^)] ^^^ds. 
Jo 
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On the other side, if <5 < 0 then €2 € R and 63 = ei, with 
3ei 7 ^ 0 , and 


2K{k) , 1 

= — ^ k = - 
2 




+ 6163 ■ 


In particular, when the nonlinear interaction is negli¬ 
gible, that is rj = 0, then the three solutions simply are 
ei = I (1 -I- p^), 62 = 63 = — I (1 -I- p^) and the period 
T actually does not depend on the initial wave function 
ipQ, but it only depends on p: 


T = 


TT 

\/l + P^ 


(7) 


In particular, in the limit of p = 0 we recover the unper¬ 
turbed beating period T = n. 

However, in general the period T depends on the initial 
wave-function, as well as on the two parameters p and rj. 
In order to estimate the dependence of the beating pe¬ 
riod T from the initial state and from the parameter p, 
corresponding to the strength of the nonlinear term, we 
consider, at Hrst, the case where the initial wave func¬ 
tion ipQ corresponds to a minimum value for the energy 
Hamiltonian R defined by ([^. As appears in Fig. [^the 
period T, corresponding to the energy minimum, actually 
depends of p for fixed value of p; for instance, the value 
of the period T at p = 0 and p = 4.2 is approximatively 
one half of the period T at p = 0 and p = 0. Only for 
large value of p we have a good agreement between the 
values of T for different values of p. 


We consider now the case where we fix the value of the 
adimensional parameters p and p and we compute the 
period T as function of the energy of the initial state. For 
argument’s sake we perform two numerical experiments 
(see Fig. [^; in the hrst one we hx p = 0.5 and p = 0.2, 
while in the second one we hx p = 2.5 and p = 4.2. For 
small values of the two parameters p = 0.5 and p = 0.2 
then the period T takes values (in the adimensional unit 
T dehned by ([^) from T = 2.754 at A = —2.117 to 
T = 2.853 dX E = 2.358, that is the period T lies in an 
interval which length is around 1.7% of the mean value 
of T, the value of T corresponding to the minimum value 
of the energy is T = 2.754. On the other hand, for larger 
values of the two parameters p = 2.5 and p = 4.2 then 
the period T takes values from T = 0.681 at E = 9.416 to 
T — 1.706 at E = —0.9586, in such a case we have that 
the period T lies in an interval which length is around 
43% of the mean value of T, the value of T corresponding 
to the minimum value of the energy is T = 1.68. Hence, 
we can conclude that for some values of the parameters 
the period T can strongly depend on the initial state. 


We close by giving the expression of the adimensional 
quantities p and p as function of physical parameters 



tl=1.2 



11 = 4.2 



FIG. 1. In this figure we plot the graphs of the period T ver¬ 
sus the adimensional parameter p for given values of p = 0.2, 
p = 1.2 and p = 4.2 and when the initial stationary state tpo 
is associated to a minimum value for the energy Hamiltonian 
R defined by ([^. Broken line correspond to the limit case 
of p = 0, in such a case the period T is given by equation 0’ 
and it does not depend on the initial wave function tpo 


in the semiclassical limit of small h; these formulas 
will be useful in order to compute the period T in a 
real device. To this end let us assume that the one¬ 
dimensional symmetric double-well potential is such that 
V{—x) = V(x) with two non-degenerate absolute minima 
points at a: = ±d, where 2d is the distance between the 
bottom of the two wells, such that 


V{±d) < 0, 


dV{±d) 

dx 


0 and 


p := 


d'^V{±d) 

dx'^ 


we assume also that the double-well potential goes to zero 
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p=0.5 ti=0.2 



p=2.5 ti=4.2 



FIG. 2. In this figure we plot the graph of the period T 
versus the values of the energy E of the initial state for given 
values of p and t]. 


for large x. The eigenvalue equation + = E(p 

admits two ground states E± = 12 w, where 17 = 
^Tiy/Ji/m [1 + 0{K)] is the ground state energy of the sin¬ 
gle trap in the limit of small with associated normal¬ 
ized eigenvectors (/ 3 ±. These eigenvectors are even and 
odd-parity functions ^p±{—x) = -Eip±{x) and the normal¬ 
ized vectors ^pr.l are given by 


fR 


V2 


and ifL 




By construction, (Pr,r{—x) = ^pl,r(x) and ipR is mostly 
localized on one of the two well (say the well with center 
at a; = +d), and ipR is mostly localized on the other well 
(say the well with center at a: = —d). In particular, 


V’Rix) = (p{x — d) and ipl{x) = (p{x + d) 

where ‘pix) corresponds to the single trap eigenvector, 
that is 

pix)=a{x;h)e-y^^"^^^ 


with 

in the semiclassical limit of small h. Hence, if we denote 
by an = [h^/mp] ^ the ground state oscillator length 


then the leading term of the parameters p and rj are given 

by 

mgd e 1 

p = - and 77 =- r= 

CO CO aijVSTT 

where co is half of the energy splitting, 2d is the distance 
between the bottoms of the two wells and e is the strength 
of the Bose-Einstein condensate. 

In conclusion; in this paper we have explored how re¬ 
late the period T of the beating motion and of the phase- 
shift of a couple of accelerated condensates separated by 
a barrier with the physical parameters of the nonlinear 
asymmetrical double-well model, mainly the strength of 
the nonlinear term, the strength of the Stark potential 
(which breaks the symmetry of the double-well poten¬ 
tial) and the splitting between the first two onsite ener¬ 
gies. We have also seen that in general the period may 
be strongly dependent of the energy of the initial state. 
Therefore, experimental determination of the physical 
constants (typically the Earth’s acceleration constant g) 
should take into account such an effect. 
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